Indication of Anisotropy in Electromagnetic Propagation over Cosmological Distances 
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We report a systematic rotation of the plane of polarization of electromagnetic radiation prop- 
agating over cosmological distances. The effect is extracted independently from Faraday rotation, 
- - - and found to be correlated with the angular positions and distances to the sources. Monte Carlo 

, analysis yields probabilistic P-values of order 10""^ for this to occur as a fluctuation. A fit yields 

■ a birefringence scale of order 10^^ ^m. Dependence on redshift z rules out a local effect. Barring 
I ' hidden systematic bias in the data, the correlation indicates a new cosmological effect. 

^ ! PACS numbers: 98.80.Es, 41.20.Jb 

' Polarized electromagnetic radiation propagating across the universe has its plane of polarization rotated by the 
' Faraday effect [Q. We report findings of an additional rotation, remaining after Faraday rotation is extracted, which 
OA i may represent evidence for cosmological anisotropy on a vast scale. 

We examined experimental data ^ on polarized radiation emitted by distant radio galaxies. The residual rotation 
^ is found to follow a dipole rule, depending on the angle 7 between the propagation wavevector k of the radiation and 
\^ a unit vector s. The rotation is linear in the distance r to the galaxy source; in sum, the rotation is proportional to 
^\ , rcos(7). This effect can not be explained by uncertainties in subtracting Faraday rotation. We focus on a statistical 
analysis of the correlation, but we have also made considerable effort to explain it in a conventional way. Unless the 
effect is due to systematic bias in the data, it seems impossible to reconcile it with conventional physics. 

Some history is useful. In 1950, Alfven and Herlofson ||l| predicted that synchrotron radiation would be emitted from 
, galaxies, with polarization perpendicular to the source magnetic field. By the mid 1960's, data began to accumulate 
on the polarization of radio waves that had traveled over cosmological distances j^]. The observables include the 
redshift z of the galaxy source, an angle ip labeling the orientation of the galaxy major axis, the percent magnitude 
of polarization p, and angles 0{X) labeling the orientation of the plane of polarization of radio waves of wavelength A. 
O ■ Experimental fits ||^ show that the angle 9{X) for a galaxy is given by 6{X) = + x- Fits to the linear dependence 
^ ' of the angle 6* on verify the presence of Faraday rotation. The fitting parameter a, called the Faraday rotation 
I measure, depends upon the magnetic field and the electron density along the line of sight Conventional Faraday 
rotation does not account for angle x, the orientation of the polarization plane after Faraday rotation is taken out; x 
is central to our analysis. 

■ On symmetry grounds, x would be expected to approximately align with the major axis angle t/i of a galaxy. This 
^ , expectation has consistently been at odds with the data. Gardner and Whiteoak j|] proposed a "two-population" 
■ - ' hypothesis, at first on the basis of 16 sources, with some sources emitting at x — "0 = 0°: and others at x — "0 = 90°- 

Clarke et al. ^ found a subset of galaxies supporting the two-population idea, but by making a severe cut consisting 
of galaxies with quite strong polarization, which eliminated most of the data {p > 11, leaving 47 of 160 galaxies.) 
The group selected as perpendicular emitters was found to be distant (high luminosity), while the parallel group was 
found to be near (low luminosity). The full data set (no strong polarization cuts) does not convincingly support 
multi-populations, and the statistical significance of conclusions is not given. The reader is warned of inconsistencies 
in the statistics: in Clarke et al. 0], the quantity A = |x — "01 is defined to be the statistic, but actually apphed when 
X — V'l < ■'''/2, while A = tt — |x — V"! is used otherwise. Carroll et al. ^ use x — "0 as a measure of residual rotation, 
a definition which neglects the other possible rotation x — "0 i (A few errors in this paper's transcription of the 
data have been corrected.) 

Birch in 1982 observed a dipole rule correlation of polarization angles and source location angles relative to an axis 
he fit from the data Birch used the acute angle between x and in a limited sample of data. The acute angle 
is an improper statistic for the observables, which are not vectors, but planes. Kendall and Young (KY) confirmed 
Birch's conclusions Q using a proper (projective) statistic. Bietenholz and Kronberg (BK) also confirmed the same 
correlation with a different analysis However, introducing a larger sample of data with sources for which redshifts 
were not known, BK then found no correlation of Birch's type. Although the subject died out after BK's negative 
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conclusion, the history is relevant inasmuch as the redshift is found to play a role. 

None of these studies addressed a correlation going like r cos(7), the cos(7) being the lowest order anisotropic effect 
that might be observed, and the factor r representing the generic dependence of birefringence on distance. The length 
r must be measured in a basic unit denoted Ag. An ansatz for the residual rotation angle (3 is 



The astronomical literature has several, sometimes inconsistent measures of angle differences. The angles x ^.nd 
ijj do not label vectors, but plane orientations; they are defined only up to multiples of tt (not 27r). Analysis should 
retain information on the sign of the differences of x and ip which probes the sense (clockwise or counterclockwise) of 
the rotation. We introduce two functions and (3^ , given by 



From Eq. the rotation /? is either positive or negative depending on the angle 7. We therefore assign a rotation 
(3 — Pf to a galaxy [i) if cos(7i) > 0, and a rotation (3 — (3~ if cos(7i) < 0. This assignment necessarily introduces 
correlations because two quadrants of the data plane of (3i and 003(7^) are excluded, a point which we discuss 
momentarily. 

By definition, (3 lies in the interval [— tt, tt]. If the residual rotation magnitude is bigger than tt for certain sources, 
then that data would be scrambled and yield low correlations. Fluctuations in the initial orientation will contribute 
noise to the analysis. Our data consists of the most complete set we have found of 160 sources in which includes 
polarization as well as position information. The radio frequency varies, typically spanning a 1-3 GHz range. Mea- 
surement uncertainties, less than 5° for ijj and typically 5° for x, were not reported on a point-by-point basis. Since 
the existing Faraday fits to Tp — aX^ + x are done using several data points, the errors on x should be much smaller 
than a few degrees. Galaxy position coordinates ||^ are given in terms of distance (r), right ascension (R.A.), and 
declination (Deck). The majority of the data comes from the northern sky, with visual magnitude of the galaxies 
between 8 and 23. We used the = 1 (critical average mass density) relation for the distance r traveled by light as 
a function of redshift z, namely r — 10^°(lightyears)[l — (1 + z)""^/^]-^, where ho — |lO~^'^(years)~ and h is the 
Hubble constant. 

To determine the validity of Eq. (Q), we computed the linear correlation coefhcient i?data for the 160 points in the 
galaxy data set of [vi cos{ji) , Pi] for trial values of the s-direction (Deck, R.A.)s sweeping out all directions. For a 
general set (xi,?/i) of N datapoints, R is defined as R ^ {Nllxiy^ - Y^XiHyi) / {[NT^xf - {T.Xif]^[NY.yf - {Y.y,f]^}. 
For each trial direction of s, we also computed 1000 correlation coefficients i?rand from 1000 copies of the original 
data set (160 x 1000 x 2 = 320,000 false data points), where the copies had rotations Pi obtained from substituting 
random ipi and Xi vnto (^. The positional part of the data, cos(7i), was not randomized. In randomizing only the 
polarizations and major axis orientations, we created a sample with the same spatial distribution of points as the data 
itself, thus allowing a scan of the polarization data without prejudice caused by the assignment, or the spatial 
non-uniformity of the data's distribution. 

For each trial s-direction, we then compared the -Rdata from the data set with the i?rand distribution, by computing 
the fraction P of computations with -Rrand ^ ^^data- In statistics, P is called a P-value; the interpretation as a 
"probability" depends on various assumptions and details of terminology. A plot of P~^ is shown in Fig. ^(a). The 
result is stable and scaled properly as we increased the number of independent trial orientations of s; the figure shows 
our finest resolution of 410 bins covering the entire celestial sphere with bins of average solid angle of about tt/IOO. 
There is a clear excess in the P~^-plot in the region s = (Deck, R.A.) = (—10° ± 20°, 20hrs ± 2hrs). 

To explore this, we cut the data to z > 0.3, roughly the most distant half of the sample (71 galaxies), to improve 
the experimental "lever arm." The correlation for the z > 0.3 set is much more dramatic; we see in Fig. 11(b) a 
well-connected cluster of more than twenty-one peaks in the region s* = (Deck, R.A.)* — (0° ± 20°, 21hrs ± 2hrs) with 
a P-value lower than we can resolve {P < 10"'^). [Several of the s-directions displayed in Fig. |l|(b) had no Monte 
Carlo events with i?rand > ^^data in the 1000-trial runs. P = 10~^ was assigned to these directions.] We call this 
procedure 1. 

For the s-direction with highest P^^-value of the full data set, the distribution of -Rrand is a Gaussian, centered 
at /i = 0.60 with a standard deviation <t = 0.032, with -Rdata = 0.66 = + 1.88cr. In contrast, in a typical "off-axis" 
direction (60°, 12hrs), the distribution is given by ^ — 0.47 and a — 0.04, with -Rdata = 0.48 = fi + 0.25(t. Returning 
to the strongly correlated data set, with z > 0.3, and for an s-direction yielding a high P~^-value, a typical /i = 0.76 
and a = 0.027, with -Rdata = 0.86 = ^J- + 3.7cr. Distributions with long tails were not seen. The spatial distribution of 
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galaxies in the sample is quite non-uniform, so that the populations assigned to /3 , and the ^ and a values, depend 
on the trial s. The P-values are therefore much more meaningful than the correlation coefficients R themselves. We 
include a few scatter plots [e.g. Figs. |l| (c), (d)], but warn that visual inspection is not very rehable. For the set of 
strongly polarized sources with polarization p> 5 (116 galaxies), the plot of versus s is almost identical to Fig. 
|l|(a). We also studied the dependence in shells of Zmm < z < Zmax containing 20 galaxies, which we considered a 
minimum number for a sensible analysis. Only in the region < z < 0.2 was P > 10~^; bin by bin, P < 10^^ in the 
region s=(0° ± 20°, 21hrs ± 2hrs) thereafter. 





FIG. 1. Inverse P-values versus anisotropy direction s, given in terms of its declination and right ascension. P is the 
fraction of galaxy sets with randomized major axis (i/'rand) and polarization (Xrand) angles that yielded a linear correlation 
coefficient -Rrand of the set [r^ cos 7^, /3i(-i/)rand, Xrand)] greater than or equal to the linear correlation coefficient iidata of the 
original set cos 7^, x)]- (a) All data (160 galaxies), (b) All data with redshift z > 0.3 (71 galaxies). The figures support 
an anisotropy direction in the region s* = (Decl., R.A.)* = (0° ± 20° , 21hrs ± 2hrs). (c) Scatter plot of /3 versus r cos 7 for the 
s-direction yielding the highest peak (arrow) of (a) for all data, (d) Scatter plot of P versus r cos 7 for the s*-direction yielding 
a typical central peak (arrow) of (b) for the cut z > 0.3. 

These studies use the correlation coefficient R appropriate to test unconstrained linear fits of the form (3 = 
iA~^rcos(7) + where A^ and S are free parameters. We found that S is consistent with zero {6 < 2°, with S 
typically « 1°) in the region of good correlation {z > 0.3, P~^ > 1000 or z > 0, P^^ > 25). We also repeated the 
entire study using R = '^Xiyi/ i^xf^yfY/''^ , which tests the hypothesis of linear correlation with the constraint 
(5 = 0, without finding significant changes. 

Another study (procedure 2) used a different order. For each random data set (again, galaxy positions were not 
randomized, as explained above), we varied s over the sphere (410 directions) to maximize i?rand- This "largest -i?iand" 
value was then recorded. A new random set was then generated, producing another "largest-i?].and-" This calculation 
was repeated more than 1000 times, to create a set of largest i?rand's. This procedure was motivated by the fact that 
there is an increased probability in procedure 1 of obtaining a fit of s to the data due to the two degrees of freedom 
of s. The crucial test is for the far-half data set with z > 0.3, which had a P-value of order 10"'^ (in procedure 1). 
For the far-half sample with z > 0.3, we found that the fraction of the largest-i?rand's that exceeded i?data(z>o.3,s=s*) 
was less than 0.006. The distribution of largest -i?iand's was characterized by a cr and /i such that i?data = M + 2.8(7. 
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In contrast, for the closest half of the data, z < 0.3, the fraction of the largest i?rand exceeding -Rdata(z<o.3,s=s*) was 
0.86, confirming that the effect "turns on" only for the most distant half of the galaxies. We found no instances of 
largest-i?rand's for the sample with z < 0.3 exceeding Rdata{z>Q.3,s=s*)- These results provide a more stringent test 
and corroborate the conclusion of procedure 1. (The authors welcome requests for additional information.) 

The average (procedure 1) best fit value is = (1.1 ± 0.08)10^^^m for a s-direction of (Decl., R.A.)* = 
(0° ± 20°,21hrs ± 2hrs) for the data with z > 0.3. For the full data set of all 160 data points, we find 
As = (0.89 ± 0.12)1025.^m for (Decl., R.A.) = (-10° ± 20°,20hrs ± 2hrs). The scale A^, approaching a billion 
parsecs, is approximately an order of magnitude larger than the largest scales observed in galaxy correlations. (Errors 
in As are the usual 1 sigma variation of uncorrelated analysis and do not refer to probabilities.) The direction s* 
appears unremarkable; although vaguely toward the galaxy center, the cluster in Fig. |l|(b) is separated from the 
galaxy center by 30°, about 30 times the apparent size of the core region of the Milky Way. 

As a consistency check, we have separately investigated whether there are strong correlations in P — A~^r, and 
P = /3ocos(7). We find nothing significant in the first case, or for the second case over the full data set. The set 
z > 0.3 does produce correlations with f3o 003(7) which we cannot readily distinguish from Eq. (|l|), since the r-values 
do not vary enough. 

As for conventional physics, the effect observed is not explained by variations on Faraday physics. While fits to 
Faraday rotation (linear in A^) represent a model and an approximation, the ratios of the radio frequency to cyclotron 
and plasma frequencies are such that the approximation is thought to have exceedingly small corrections. Observers 
also make corrections for systematic errors and take into account the effects of the Earth's ionosphere. Consulting 
with the original observers Q does not yield any suggestions for bias that would imitate the signal we observe. 
We have nevertheless questioned whether a local effect of the galaxy, via some unanticipated conventional physics, 
might account for our correlation. The fact that the correlation is seen for z > 0.3, but not z < 0.3, rules out a 
local effect. (Several observational groups contributed to the z > 0.3 data, closing the loophole that one particular 
analysis might contain bias, as best we can determine.) Again, strong fields at the source might generate unexpected 
initial polarization orientations, or upset the Faraday-based fits, and this could plausibly depend on z. But since 
the correlation is observed in cos(7), any population-based explanation requires an unnatural, if not impossible, 
conspiracy between distant sources at widely separated zenith angles. One is left, then, with the option of unknown 
systematic bias for the large-z set, or accepting the possibility that the correlation is a real physical effect. For the 
latter, one must arbitrarily invoke coherence on outrageously vast distances, perhaps organized by electromagnetic or 
other interactions in the early universe, or contemplate new physics. 

If we take the data at face value as indicating a fundamental feature of electrodynamics, gauge invariance severely 
limits the possible couplings of the vector potential A^^ and the electromagnetic field strength tensor F^" = d'^A" — 
d'^A'^ to any background vector s'^. The unique derivative expansion (units are h — c — 1) for terms in the effective 
action Sgff is 

Seff = J d'xi^-^F^^F^^'^ + ^A^'e^^'-'^^F^^A^spy (3) 

suppressing higher derivative terms which would contribute to short-distance effects. The dispersion relation for this 
theory at lowest order in A~^ is k± = to ± ^Aj^ 003(7), where = (0, s) in the coordinate system where w and k are 
measured. Rotation of the plane of polarization comes from differences in propagation speed between the two modes; 
the difference ^(fc+ — fc_) = d(3/dr is a measure of the polarization plane rotation (3 per unit path length r, yielding 
a rotation (3 coinciding with Eq. (p. 

The interpretation of the parameters s'' depends on their physical origin and their transformation properties. If s'^ 
is odd under time reversal, and its space part is a pseudovector under parity, then Eq. (|^) preserves these symmetries 
separately. If s'^ is also spacelike, as we have assumed, one might associate the vector swith an intrinsic "spin axis" of 
an anisotropic universe. Terms of order are dropped consistently if our reference frame coincides well enough with 
a rest frame of s''; s° can also be added as a separate parameter. Because the new terms have one less derivative than 
the standard ones, these terms have no effect (by power counting) on high energy questions such as renormalizability. 
The scale As represents a fundamental length scale in the modified electrodynamics due to Eq. (||), and should not 
be confused with a "photon mass" , which violates gauge invariance. By converting the length A^ to a mass scale 
lj.c^ = A~^hc, a value of n — 10~^^eV/c^ is found, which is 10^^ times smaller than the photon mass limit of Chibisov, 
and lO^'' times smaller than Goldhaber and Nieto's 

Ni 1^ obtained (^) from covariance arguments; independently it was found from quantum adiabatic arguments 
(Ralston §). The latter was our initial motivation, used to predict the correlation (|^), which led to this investigation. 
The curious history came to light much later. If is treated as dependent on x, then it must be a gradient: 
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Aj^s'^ = ^d'^cj), where g is a coupling. The theory is then related to axions or similar pseudoscalar fields (p with 
coupling g. If a new field is proposed, there should be further observational consequences. Our study of a spacelike 
is unconventional compared to other work , but there is every reason to think that the observed correlation can be 
consistent with axion-type domain walls 0, or other condensate structures. More data on the observables is needed, 
especially from the southern sky. The crucial issues of conventional explanations and experimental systematics merit 
scrutiny from a broad community. From a scientific standpoint, we report what we find, given the data that exist. We 
find that the data contain a correlation indicating cosmological anisotropy in electromagnetic propagation. Further 
study may be able to determine whether or its counterparts invoking new fields might be a valid description of 
electromagnetism on the largest scale. 

B. Anthony-Twarog, K. Ashman, C. Bird, H. Rubinstein, B. Cox, and an anonymous referee made helpful sugges- 
tions. DOE Grant Number DE-FG02-85ER40214, NSF Grant Number PHY94-15583, and the K*STAR program 
provided support. 
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